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LONG-TERM  GOALS 

The  long-term  goals  are  to  contribute  improvements  in  current  physical  understanding  and  modeling  of 
interfacial  processes  fundamental  to  air-sea  interaction  fluxes,  particularly  those  involving  wave 
breaking  and  spray  droplet  production.  These  advances  will  improve  the  reliability  of  operational  sea 
state  and  ocean  weather  forecasting  models,  particularly  for  severe  to  extreme  sea  states. 

OBJECTIVES 

This  project  seeks  to  improve  the  reliability  of  air-sea  interfacial  flux  parameterizations  in  coupled  sea 
state/marine  weather  forecasting  models,  with  a  particular  focus  on  refining  and  incorporating  the  role 
of  wave  breaking  and  sea  spray  in  severe  conditions.  The  approach  adopted  is  to  develop  more  realistic 
parameterizations  for  breaking  occurrence  and  strength,  and  sea  spray/spume  source  functions  and 
validating  them  in  test-bed  models.  The  end  goal  is  implementing  these  improvements  in  a  coupled 
COAMPS/WaveWatch  III  model  for  operational  use. 

APPROACH 

Our  approach  is  to  build  substantially  on  our  accumulated  expertise  in  sea  surface  processes  and  air-sea 
interaction  (Banner)  and  numerical  weather  modelling  (Leslie)  to  identify  and  close  fundamental 
knowledge  gaps  in  order  to  improve  the  modeling  accuracy  for  severe  marine  meteorological  events 
such  as  hurricanes.  Such  improvements  depend  critically  on  gaining  a  more  complete  understanding  of 
severe  sea  state  phenomena  linked  to  wave  breaking,  due  to  their  increased  surface  drag  and  allied  air- 
sea  flux  enhancements,  including  sea  spray  production. 

Quantifying  the  distribution  of  wave  breaking  events  has  progressed  substantially  within  this  project  as 
a  result  of  an  intensive  collaboration  with  D.  Farmer  and  J.  Gemmrich  (IOS,  Canada)  in  FY01 .  From 
their  storm  wave  datasets,  a  robust  threshold  behavior  was  identified  for  wave  breaking  at  different 
scales  in  terms  of  the  wave  spectral  saturation  level  B(/)  (Banner  et  al,  2002,  hereafter  BGF02).  Here 
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B (/)  is  given  by  B(/)  =  (2;z)4/5F(/)/2g2,  where  F//)  is  the  wave  energy  spectrum.  This  result  has  proven 
very  useful  in  our  ongoing  effort  to  parameterize  wave  breaking  spectrally  in  wind-wave  models,  and 
underpins  our  new  spectral  source  term  formulations  for  the  wave  energy  dissipation  rate.  It  is  also 
being  used  in  the  spray/spume  production  term  formulated  in  collaboration  with  C.  Fairall  (NOAA). 

Model  development  in  this  project  has  focused  on  refining  full-bandwidth  computations  of  the 
directional  wave  spectrum  and  its  tail  region  using  ‘exact’  versions  of  the  nonlinear  wave-wave 
interaction  source  term  in  the  radiative  transfer  equation  for  the  wave  field.  We  have  concentrated  on 
fetch-limited  and  duration-limited  growth  cases,  since  the  bulk  of  observational  data  exists  for  these 
cases.  This  ensures  that  modeled  spectral  saturation  levels  (and  dissipation  rates)  are  consistent  with 
observed  levels,  which  is  essential  for  predicting  spectral  breaking  wave  properties  and  reliable 
calculation  of  the  enhanced  air-sea  fluxes  associated  with  wave  breaking.  This  is  needed  to  underpin 
future  operational  versions  of  coupled  atmosphere -wave-ocean  models. 

We  have  focused  strongly  on  formulating,  implementing  and  refining  modeling  strategies  for:  (a) 
extracting  the  relevant  wave  breaking  parameters  and  (b)  calculating  wind  stress/roughness  length 
enhancements  and  updating  the  surface  layer  winds  accordingly.  For  (a),  we  have  significantly  refined 
our  previously  reported  capabilities  for  calculating  the  spectral  density  of  mean  breaking  wave  crest 
length/unit  area.  This  is  a  primary  goal  of  this  project,  as  this  quantity  is  central  to  the  prediction  of 
breaking  wave  enhancements  to  the  wind  stress,  and  development  of  a  spray/spume  source  function 
based  on  sea  state  rather  than  the  wind  strength.  Progress  has  been  made  by  refining  the  form  of  Sds 
proposed  by  Alves  and  Banner  (2002)  [hereafter  AB02]  and  on  advances  in  predicting  wave  breaking 
probabilities  at  different  wave  scales  in  BGF02.  In  particular,  the  AB02  form  of  Sds  was  upgraded  to 
incorporate  the  observed  BGF02  breaking  saturation  threshold.  Also,  various  refinements  were 
introduced  to  its  spectral  distribution  in  order  to  provide  a  much  better  match  to  the  wind  input  source 
function  Sin  at  higher  wavenumbers,  and  to  recently  published  spectral  wave  breaking  observations  of 
the  spectral  density  of  breaking  crest  length/unit  area.  To  progress,  we  have  had  to  invest  very 
considerable  effort  evaluating  various  versions  of  the  ‘exact’  Sni  code,  wind  input  source  tenns  Sin,  and 
propagation  schemes  in  the  radiative  transfer  evolution  equation,  in  order  to  ensure  accuracy  and 
minimise  computational  instabilities  that  can  develop  at  higher  wavenumbers. 

WORK  COMPLETED 

1.  Wind- wave  evolution  modeling. 

A.  uncoupled  study 

During  FY05  we  continued  our  main  effort  on  hurricane  wind-wave  model  development,  with  explicit 
computation  of  wave  breaking  properties.  Our  capability  for  computing  the  directional  wave  spectrum 
was  refined  further,  with  particular  emphasis  on  the  directional  spectral  tail  region  out  to  high 
frequencies.  Attention  was  focused  on  extending  the  wind  speed  from  the  commonly  observed  7-14 
m/s  wind  speed  range  out  to  30  m/s  and  beyond.  Direct  comparisons  of  model  results  with 
observations  were  made  against  ID  transect  wave  spectra  and  the  breaking  crest  length  spectra  of 
Melville  and  Matusov  (2002,  hereafter  MM02)  and  Gemmrich  (2005),  as  well  as  the  directional 
spreading  data  of  Hwang  et  al.  (2000).  For  the  very  high  wind  speeds  representative  of  hurricane 
winds,  our  aim  was  to  achieve  stable  model  behavior,  especially  in  the  spectral  tail  region. 
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In  the  ‘uncoupled’  modeling  phase,  the  wind  input  source  term  S,n  was  based  on  a  modified 
parameterisation  of  Yan  (1987),  which  embodies  the  traditionally-observed  growth  rate  dependences 
on  wind  parameters  linear  near  the  spectral  peak  and  quadratic  in  the  spectral  tail.  It  also  provides 
energy  input  to  the  wind  from  waves  outrunning  the  wind.  These  computations  also  used  a  state-of-the- 
art  version  of  the  ‘exact’  nonlinear  source  term  Sni  (D.  Resio,  private  communication).  Our  spectral 
wave  energy  dissipation  rate  source  term,  5*,  has  a  primary  dependence  on  wave  nonlinearity  through 
its  functional  dependences  on  the  spectral  saturation.  The  uncoupled  model  performance  was 
investigated  by  running  the  spectral  wave  evolution  equation  for  duration  and  fetch-limited  conditions, 
computing  the  spectrum  over  a  wide  bandwidth,  typically  out  to  several  Hz. 

Throughout  these  calculations,  the  results  for  integral  wave  spectral  properties  such  as  mean  energy 
and  peak  frequency  achieved  a  close  correspondence  with  the  observed  evolution  behavior  in  Kahma 
and  Calkoen  (1992),  as  shown  in  previous  annual  reports.  We  also  obtained  very  close  model 
agreement,  both  level  and  spectral  slope,  with  the  high  wavenumber  tail  of  the  1-D  transect  spectrum 
observed  by  MM02,  also  shown  in  previous  annual  reports.  Lastly,  modeled  directional  spreading 
properties  compared  very  favorably  with  the  data  of  Hwang  et  al.  (2000)  and  the  more  recent  high 
wind  speed  results  from  the  recent  GOTEX  experiment  (Melville/Friehe),  in  regard  to  both  the  bimodal 
shape  at  higher  wavenumbers  and  the  mean  spreading  width  against  wavenumber. 

On  this  basis,  we  proceeded  to  make  extensive  calculations  of  the  spectral  density  of  breaking  crest 
length  per  unit  area,  and  how  these  breaking  wave  spectra  change  with  wind  speed  and  wave  age. 

While  awaiting  CBLAST  hurricane  data  from  the  SIO  camera  deployments,  the  field  observational 
results  of  MM02  and  Gemmrich  (2005)  in  the  7-14  m/s  wind  speed  regime  have  been  used  to  validate 
our  model  results. 

From  the  spectral  evolution  calculations,  of  particular  importance  is  the  spectral  distribution  of  the 
dissipation  rate  S*.  The  precise  link  between  S*  and  the  spectral  density  of  breaking  crest  length/unit 
area  A(c)  is  a  current  research  topic  of  considerable  complexity. 

The  assumption  that  wave  breaking  dissipation  is  localized  spectrally  is  not  well-established,  i.e.  that 
spectral  dissipation  at  scale  c  (or  k )  is  only  associated  with  breaking  waves  of  that  scale.  If  localized 
spectral  dissipation  is  assumed,  following  Phillips  (1985),  calculation  of  A(c)  follows  from: 

SdS(c)  =  b(c5/s)  A(c)  (1) 

The  non-dimensional  coefficient  b  reflects  the  strength  of  breaking,  but  its  spectral  dependence  is  not 
known  (Melville,  1996)  and  needs  to  be  assessed  observationally.  As  there  are  no  available  direct 
measurements  of  spectral  dissipation  rate  5*,  the  available  field  measurements  of  A(c)  of  MM02  and 
Gemmrich  (2005)  can  be  used  to  investigate  whether  our  S*  model  provides  a  self-consistent  and 
robust  fonnulation  under  the  localized  spectral  dissipation  assumption. 

There  are,  however,  several  known  non-local  sources  of  dissipation,  such  as  attenuation  of  short  waves 
both  by  longer  breaking  waves  and  also  through  strong  nonlinear  interaction  with  non-breaking,  steep 
longer  waves.  In  strongly  forced  sea  states  typical  of  hurricanes,  these  spectrally  non-local  effects  are 
likely  to  be  significant.  Our  model  development  effort  has  been  extended  to  include  such  effects. 
Recent  progress  is  described  below  under  ‘New  Results’. 
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B.  coupling  wind  profile  and  sea  surface  aerodynamic  roughness  changes 

It  is  well-recognized  that  forecasting  wind-waves  involves  strong  coupling  between  the  wind  and  wave 
fields,  and  that  the  growth  of  waves  is  accompanied  by  corresponding  changes  in  the  marine 
atmospheric  boundary  layer.  Embodying  the  complex  physics  of  this  coupling  is  an  essential  part  of 
progress  in  improved  forecasting.  To  address  this  challenge,  we  have  developed  a  Fortran  code  that 
allows  the  iterative  computation  of  the  wind  stress  and  wave  spectrum  in  response  to  changes  due  to 
wave  drag,  including  breaking-induced  contributions.  Our  approach  is  based  on  the  calculation  of  wave 
breaking  properties  from  the  wave  spectrum,  the  consequent  enhanced  wave  drag  and  the  adjustment 
of  the  aerodynamic  roughness  length  zo  in  the  assumed  logarithmic  mean  velocity  profile  for  the 
surface  layer  wind  field.  In  turn,  the  wave  spectrum  is  modified  in  response  to  the  updated  wind 
profile.  This  has  been  incorporated  into  our  coupled  test-bed  model,  and  we  have  been  investigating  its 
performance  in  great  detail  over  a  wide  range  of  wind  speeds  up  to  hurricane  strength.  An  important 
underlying  quantity  is  A(c),  and  the  available  data  is  crucial  in  validating  the  breaking  wave  extraction 
algorithm  in  the  model.  Representative  coupled  model  results  are  described  below  under  ‘Results’. 

2.  Sea-state  dependent  spume/spray  droplet  source  term 

In  collaboration  with  C.  Fairall  and  W.  Asher,  an  exploratory  laboratory  experiment  was  staged  in 
2003.  The  aim  was  to  seek  a  dependence  of  spume  droplet  flux  on  breaking  wave  properties  rather  on 
the  wind  speed.  This  approach  underpins  a  new  model  of  spume  droplet  production  in  terms  of  the 
surface  dissipation  rate  of  breaking  waves  developed  with  C.  Fairall.  Droplet  measurements  were  made 
in  wind  velocities  in  the  range  Uio~  20-35  m/s,  with  paddle-initiated  wind  waves  with  frequencies  of 
1.36  -1.6  Hz.  Salinity  was  varied  from  0  to  25  ppt.  The  surface  properties  were  measured  with  a  small 
directional  array  of  wire  impedance  gages  each  sampled  at  over  2000  Hz.  An  analysis  was  undertaken 
to  extract  a  consistent  measure  of ‘surface  disturbance  energy  level’.  The  results  suggest  a  fundamental 
link  between  the  spray/spume  flux  and  the  strength  of  breaking  events.  These  findings  indicate 
extrapolation  to  oceanic  scales  is  a  potentially  fruitful  future  challenge.  A  parametric  model  based  on 
this  approach  has  been  developed  for  testing  in  full  physics  models.  A  manuscript  describing  both  the 
observations  and  the  model  is  nearing  completion.  More  details  were  given  in  the  FY04  report. 

NEW  RESULTS 

Coupled  model  results  for  hurricane  wind  speeds  (30  m/s) 

Previous  annual  reports  for  this  project  show  the  excellent  correspondence  between  modeled  and 
observed  evolution  for  the  key  integral  and  spectral  tail  region  properties,  based  on  our  generic 
saturation-based  fonn  for  Sds  and  the  wind  input  source  term  Sin  due  to  Yan  (1987).  The  exact  form  of 
the  nonlinear  spectral  transfer  term  S„i.  These  results  were  obtained  at  lower  wind  speeds  (7-15  m/s) 
assuming  no  coupling  with  the  overlying  wind  profile  and  were  described  in  previous  annual  reports. 
Different  wind  input  functions  can  be  accommodated  by  adjusting  the  exponents  of  the  saturation- 
dependent  tenns  in  Sds- 

In  the  uncoupled  model  phase,  we  also  carried  out  complementary  computations  of  A(c),  the  spectral 
density  of  breaking  crest  length/unit  area  of  sea  surface,  expressed  as  a  function  of  the  wave  phase 
speed  c.  The  spectral  fonn  of  A  was  based  on  the  local  spectral  breaking  assumption  using  eqn  (1) 
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above,  and  our  calculations  indicated  that  the  main  balance  in  the  tail  is  between  Sjn  and  Sds,  with  Sni 
secondary.  The  wind  speeds  investigated  were  those  (7.2  and  13.6  m/s)  of  the  old  wind  sea 
observations  of  MM02.  When  scaled  by  the  cube  of  the  friction  velocity,  the  model  results  reproduced 
well  the  close  to  cubic  wind  speed  dependence  observed  by  MM02,  and  the  c'6  spectral  dependence 
obtained  with  a  constant  value  of  the  breaking  strength  coefficient  b  in  (1)  matched  the  theoretical 
predictions  of  Phillips  (1985)  f  or  A(c)  in  the  equilibrium  subrange  towards  the  spectral  peak. 

For  the  coupled  model,  we  used  the  latest  available  form  of  Sin  from  the  Lake  George  study  (Donelan 
et  al.,  2005).  The  model  Sm  is  fonnulated  using  the  wind  at  a  height  of  one-tenth  of  the  local  wave 
length,  U(k!\  0),  rather  than  u*.  and  includes  an  explicit  computation,  and  updating,  of  the  aerodynamic 
roughness  length  due  to  the  enhanced  wave-coherent  component  of  the  wind  stress,  as  described  in 
section  IB  above. 

We  have  also  generalized  the  localized  spectral  breaking  assumption  used  in  our  uncoupled  study,  and 
now  seek  to  incorporate  parametrically  the  dissipation  due  to  additional  known  mechanisms.  These 
include  the  removal  of  short  wave  energy  by  the  cumulative  effect  of  longer  breaking  waves  (e.g. 
Banner  et  al,  1989),  and  also  by  strong  nonlinear  interaction  with  steep  dominant  waves  (e.g.  Phillips 
and  Banner,  1974).  The  spectral  distribution  A(c)  is  now  calculated  from  the  residual  local  dissipation 
due  to  breaking  at  scale  c,  after  the  non-local  contributions  to  dissipation  at  that  spectral  scale  have 
been  subtracted.  The  following  figures  show  the  key  data  products  from  a  duration-limited  run  where 
the  wind  speed  at  200m  was  set  to  35  m/s.  As  the  spectrum  evolves,  the  model  calculates  the  changing 
surface  layer  winds,  in  response  to  the  sea  state,  which  in  turn  modifies  the  near-surface  wind  profile. 

The  following  figures  show  key  aspects  of  the  wave  field  during  the  evolution.  These  include  key  wave 
breaking  spectral  properties,  in  addition  to  the  familiar  wave  spectra.  The  final  figure  shows  the 
evolution  of  the  wind  stress  constituents  -  non-breaking  wave  form  drag,  additional  breaking  wave 
form  drag,  and  viscous  frictional  drag. 
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Figure  1.  Evolution  of  the  non-dimensional  wave  energy  (upper  panel)  and  non-dimensional 
dominant  wave  frequency  (lower  panel)  against  non-dimensional  time.  The  background  dashed 
curve  is  the  trend  of  available  field  observations,  gathered  for  lower  wind  speeds. 
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Figure  2.  One-dimensional  transect  spectrum  in  the  dominant  wave  (k() 
direction  for  mature  seas,  showing  a  k{3  behavior.  The  background  curve  was 
measured  by  MM02  for  7-13  m/s  winds.  The  present  tail  spectral  level,  close  to 
3  times  higher,  indicates  a  near-  linear  dependence  on  wind  speed. 
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Figure  3.  (a)  azimuth-integrated  spectral  saturation  against  distance  from  the  spectral  peak 
wavenumber,  for  inverse  wave  age  Uk/cp  ~1.1  (b)  corresponding  normalized  saturation  formed 
from  (a)  by  division  by  the  mean  spreading  width  in  radians  (c)  directional  distribution  of  energy  at 
different  wavenumbers  relative  to  the  spectral  peak  (d)  variation  with  wavenumber  of  the  mean 
spreading  width  in  degrees,  compared  with  observations  reported  by  Hwang  et  al.  (2000).  The 
calculations  match  the  observed  peaked  off-wind  nature  and  strong  broadening  of  the  direction 
spreading  towards  the  shorter  wave  scales  reported  in  recent  field  studies. 


Figure  4.  Evolution  of  the  normalized  spectral  density  of  breaking  crest  length/unit  area  (A) 
against  the  wave  speed.  The  legend  shows  the  inverse  wave  age  UifCp.  The  normalization  is  based 

onU  w ,  following  MM02,  and  the  dashed  curve  shows  their  data  for  lower  wind  speeds  and  fully 

developed  seas.  The  enhanced  A  levels  at  the  evolving  spectral  peak  seen  in  these  younger  wind  sea 
calculations  tend  to  disappear  towards  full  development. 
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Figure  5.  Evolution  of  the  modeled  distribution  of  breaking  probability  with  wave  frequency,  based 
on  a  half-power  extrapolation  of  the  open  ocean  storm  results  of  Banner  et  al.  (2002)  for  winds  up 
to  gale  force.  The  legend  shows  the  inverse  wave  age  Un/Cp.  Note  the  strong  change  in  breaking 
probability  as  the  sea  ages  at  these  high  wind  speeds,  with  significant  levels  of  breaking  probability 

forecast  even  when  the  seas  are  mature. 
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Figure  6.  Evolution  with  wave  age  of  the  normalized  wind  stress  (the  sea  surface  drag  coefficient). 

It  decreases  from  about  0.0037  to  0.0033  as  the  wind  sea  ages.  At  full  development,  the  projected 
value  is  around  0. 0031.  The  normalized  components  of  the  wind  stress  are  also  shown,  with  the  non¬ 
breaking  wave  stress  increasing,  while  the  breaking  wave-induced  stress  decreases.  The  viscous 
stress  level  is  insignificant  for  these  near-hurricane  winds. 
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Overall,  these  model  results  for  near-hurricane  wind  speeds  are  very  encouraging,  but  at  this  stage  it  is 
essential  to  perform  model  validations  based  on  comparison  with  the  CBLAST  Hurricane  project  and 
other  severe  wind  observation  data  sets,  which  are  presently  in  progress.  In  the  interim,  further 
refinements  of  Sds  have  been  fonnulated  and  are  being  investigated. 

IMPACT/APPLICATIONS 

Enhanced  scientific  understanding  of  severe  sea  state  air-sea  interfacial  processes,  particularly  wave 
breaking  and  spray/droplet  production  rates,  will  provide  more  reliable  parameterizations  of  these 
processes  and  closely  related  air-sea  fluxes  when  introduced  into  operational  forecast  models.  These 
improved  parameterizations  will  increase  the  reliability  of  operational  sea  state  and  marine 
meteorological  forecasts,  especially  during  severe  marine  weather  conditions.  Of  particular  benefit  will 
be  the  capability  to  provide  routine  forecasts  of  occurrence  rate  of  dangerous  breaking  waves. 

RELATED  PROJECTS 

The  ONR  project  Source  Term  Balance  for  Finite  Depth  Wind  Waves  (Young,  Banner  and  Donelan) 
includes  a  strong  focus  on  the  influence  of  steep  waves  and  wave  breaking  on  the  wind  input  source 
function  in  strongly  forced,  constant  depth,  shallow  water  environments.  This  data  set  has  been 
analyzed  and  initial  papers  on  the  methodology  and  results  have  appeared  in  refereed  journals,  with 
other  papers  in  preparation  [e.g.  Donelan  et  al,  2004;  Banner  et  al.,  2005].  In  the  present  coupled 
modeling  effort,  we  have  implemented  the  fonn  of  wind  input  source  function  Sin  from  that  project. 
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